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ABSTRACT 

Despite its potential reactivity due to ring strain, ethylene oxide (C-C2H4O) is a com- 
plex molecule that seems to be stable under the physical conditions of an interstellar 
dense core; indeed it has been detected towards several high-mass star forming re- 
gions with a column density of the order of 10 13 cm -2 (Ikeda et al. 2001). To date, 
its observational abundances cannot be reproduced by chemical models and this may 
be due to the significant contribution played by its chemistry on grain surfaces. Re- 
cently, Ward & Price (2011) have performed experiments in order to investigate the 
surface formation of ethylene oxide starting with oxygen atoms and ethylene ice as 
reactants. We present a chemical model which includes the most recent experimen- 
tal results from Ward & Price (2011) on the formation of C-C2H4O. We study the 
influence of the physical parameters of dense cores on the abundances of C-C2H4O. 
We verify that ethylene oxide can indeed be formed during the cold phase (when the 
ISM dense cores are formed), via addition of an oxygen atom across the C=C double 
bond of the ethylene molecule, and released by thermal desorption during the hot core 
phase. A qualitative comparison between our theoretical results and those from the 
observations shows that we are able to reproduce the abundances of ethylene oxide 
towards high-mass star-forming regions. 

Key words: astrochemistry - stars:formation - ISM:abundances - ISM: molecules. 



1 INTRODUCTION 

Ethylene oxide (C-C2H4O) is the simplest epoxide, a ring- 
shaped organic compound involving an oxygen atom bonded 
to two carbon atoms. Since it is believed to lead to the syn- 
thesis of amino acids (Cleaves 2003) and hence potentially 
influence early pre-biotic pathways (Miller & Schlesinger 
1993), C-C2H4O may play an important role in the inter- 
stellar medium (ISM). 

First detected towards the Sgr B2N molecular cloud as- 
sociated with 10 rotational transitions (Dickens et al. 1997), 
this molecule has been investigated in a more extended 
frequency range (662 lines between 262-358 GHz) (Pan et 
al. 1998) and more recently Bernstein & Lynch (2009) at- 
tributed the unidentified IR band features discovered by 
Gillett et al. (1973) to ethylene oxide and to another 
small carbonaceous molecule, cyclopropenylidene. Ikeda et 
al. (2001) estimated the C-C2H4O fractional abundances by 
observing 20 massive star-forming regions and two dark 
clouds. In fact, ethylene oxide was detected towards 10 of 
these sources but it was not found in TMC-1 and TMC- 
RNH3) dark clouds, in contrast to one of its isomers, ac- 
etaldehyde (CH3CHO) which seems to be ubiquitous. Ikeda 
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et al. (2001) also investigated the relationship between c- 
C2H4O abundances and the dust temperature, but without 
finding any correlation. On the other hand, the evidence of 
similar trends in the column densities of all C2H4O isomers 
and ethanol identified the latter as a possible precursor of 
ethylene oxide and its isomers. Indeed, a gas-phase pathway 
for ethylene oxide formation has been proposed (Dickens et 
al. 1997) with ethanol as one of the reactants: 

CH 3 + + C 2 H 5 OH -)■ C 2 H 5 + + CH 4 

C 2 H 5 + + e~ -> C 2 H 4 + H 

where the formula C2H 4 does not distinguish among 
the different isomers. However, Ikeda et al. (2001) performed 
simulations of c-C2H 4 chemical evolution by including only 
pure gas-phase reactions at different temperatures and vol- 
ume densities and they found that the calculated abun- 
dances of ethylene oxide decreased when the gas-kinetic tem- 
perature was raised. The results of Ikeda et al. (2001) sug- 
gested that other reactions such as grain surface reactions 
may be important for the synthesis of ethylene oxide. 

Another gas-phase pathway for the formation of ethy- 
lene oxide has been proposed by Turner & Apponi (2001) 
which studied the reaction chain to produce C-C2H4O and 
one of its isomers, vinyl alcohol (CH2CHOH), starting with 
the interaction between CH3 + and H2CO; unfortunately, 
none of the intermediates of the latter reaction have been 
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identified yet in the ISM. The abundance ratios among the 
three isomers, acetaldehyde, ethylene oxide and vinyl alco- 
hol, seem to be very uncertain (Bennet et al. 2005), but 
generally CH3CHO was found to be the most abundant 
and CH2CHOH has been predicted to have a lower abun- 
dance than ethylene oxide in the experiments by Bennet et 
al. (2005). This comparison could be explained in the light 
of the fact that the formation of ethylene oxide occurs via 
direct addition of the oxygen across the C-C double bond 
of ethylene, while vinyl alcohol requires an insertion of an 
oxygen into the C-H bond. Hudson & Moore (2003) sug- 
gested a solid state route for the formation of CH2CHOH 
involving exposure of H20:C2H2 ices to protons and pho- 
tons, but the first proposed surface scheme for ethylene ox- 
ide formation dates back to Charnley (2004) who identified 
acetaldehyde (CH3CHO) as a possible precursor. The con- 
tribution of grain-surface reactions to the production of the 
three isomers with formula C2H4O has also been investi- 
gated by Bennet et al. (2005): they presented the results of 
electron irradiation on a mix of carbon dioxide and ethylene 
ice, under interstellar and cometary conditions, to highlight 
the effect of non-thermal processing of ices. 

Recently, an alternative solid state pathway for the for- 
mation of ethylene oxide has been investigated by Ward 
& Price (2011) who studied the surface reactivity between 
C2H4 and thermal O atoms under astrophysically relevant 
conditions. They co-deposited the reactants on a graphite 
surface at various temperatures in the 12-90 K range. In or- 
der to identify the products formed, they recorded desorp- 
tion spectra and they also selectively ionized the three differ- 
ent isomers with a tunable dye laser. Finally, they performed 
a similar laboratory study with C3H6 as reactant, but unfor- 
tunately although C3H6O was observed they were unable to 
identify the isomeric products of the reaction due to the large 
number of possible molecules with formula C3H6O. Despite 
these limitations, the authors assumed that the majority of 
the C3H6O formed is propylene oxide by following the same 
mechanism as for the case of ethylene oxide production. Un- 
like C-C2H4O, C-C3H6O has not yet been detected in the 
ISM since it is believed to isomerize to propanal or acetone, 
which have both been observed in the ISM (Hollis et al. 2004 
and Combes et al. 1987, Snyder et al. 2002). To date, the 
observed abundances of ethylene oxide have not yet been 
reproduced by chemical models; in fact, to our knowledge, 
none of the grain-surface pathways for formation for ethy- 
lene oxide have been included in the models that simulate 
the star formation process. 

In this paper, we update the UCL_CHEM gas-grain 
chemical model with the reaction rates for the production of 
this species as evaluated by Ward & Price (2011). By follow- 
ing the same logic as these authors, we extended our study 
to include the formation of propylene oxide in order to ver- 
ify if C3H6O can be produced under ISM conditions even 
without having any observational evidence of its presence 
in the interstellar space. The paper is organized as follows: 
Section 2 gives details of the model; Section 3 considers the 
astronomical implications of our study; Section 4 contains a 
discussion of our results and it also shows a comparison be- 
tween our theoretical column densities and those taken from 
the observations. Finally, we give our conclusions in Section 
5. 



Table 1. Surface reactions investigated by Ward & Price (2011) 





"Reaction 


a (cm 2 mol 1 s 


r ) at 20K 


mO - 


h mC 2 H 4 mc-C 2 H 4 


1.20x10" 


-19 


mO - 


h mC 3 H 6 -> mc-C 3 H 6 


3.40x10" 


-19 



a The m before the molecular formulas stands for mantle 
2 THE CHEMICAL MODEL 

UCL_CHEM is a time-depth dependent gas-grain chemi- 
cal model which simulates the abundances of molecules in 
the ISM and star-forming regions. The first version of the 
code was developed by Viti & Williams (1999) and revis- 
ited by the same authors in 2004 (Viti et al. 2004). Our 
model has been extensively described in the literature (see 
Occhiogrosso et al. 2011 for more details) so here we only 
highlight the main points. 

The process of star formation is simulated in two steps: 
during Phase I the gas collapses to a dense core from a dif- 
fuse and atomic state; at this stage the gas freezes onto the 
grains and surface reactions occur. During Phase II, the gas 
and the dust surrounding the central star heat up, causing 
the sublimation of the icy mantles. The thermal desorption 
can be set up to occur instantaneously or as a function of 
time as in Viti et al. (2004). In the latter case, molecules 
are classified in five different categories as described by the 
experimental study of Collings et al. (2004) : (i) CO-like; (ii) 
H20-like; (iii) intermediate; (iv) reactive and (v) refractory. 
In the present work, we consider C-C2H4O as HbOTike (it 
will co-desorb when the H2O ice desorbs) based on the fact 
that in the Ward & Price (2011) experiments, its desorp- 
tion profile shows only one peak (as for the water desorp- 
tion profile) and its dipole moment is big enough to assume 
that the species will have a lower desorption energy from 
the surface (in that specific case, graphite) than from it- 
self. In our model the temperature of the gas is derived as a 
function of the luminosity (and therefore the age) of the ac- 
creting protostar through a simple power law (see Viti et al. 
2004 for details). The model includes two networks: the gas- 
phase network is based on the UMIST database (Woodall et 
al. 2000) (even if some coefficients have been modified with 
those from the most recent database KIDA (Wakelam et al. 
2009)) with some updates discussed later (see Section 3), 
and a grain-surface network with reactions whose rates are 
experimentally determined where possible. 

For this study, we include in our surface reaction net- 
work the formation of C-C2H4O and C-C3H6O on the grain 
surface as investigated by Ward & Price (2011) (see Table 
1). In order to analyze their data, these authors developed a 
simple kinetic model in which molecules can react via Eley- 
Rideal (ER) or Langmuir-Hinshelwood (LH) mechanisms. 
In the latter case, both reactants adsorb and thermalize on 
the surface before a diffusion reaction occurs. We consider 
the reactions to occur via the LH mechanism that seems 
of greater astrophysical relevance as stated by Awad et al. 
(2005) and by the authors themselves (see their paper for 
the details). 

Since our code accounts for reactions that occur in 
three dimensions (i.e. gas-phase reactions), we need to trans- 
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Table 2. List of test chemical models. 



Test 


mH 2 (%) 


mO (%) 


mOH (%) 


Testl 


99.5 


0.5 





Test2 


90 


1 


9 


Test3 


90 


2 


8 


Test4 


90 


5 


5 


Test5 


90 


10 





Test6 


80 


20 





Test7 


70 


15 


15 



form the surface rate constants of Ward and Price discussed 
above. Therefore we first evaluate two quantities: the vol- 
ume of ice per grain and the number density of grains. In 
the latter case we assume a dust/mass ratio of 0.01 and a 
dust density of 2.5 g cm" 3 . Thus, we obtain a grain number 
density of 3.8xl0~ 5 cm" 3 . In order to estimate the volume 
of ice we consider the grains as spheres and we assume a 
thick ice layer of 0.3 urn (Collings et al. 2003). Since the 
major partner in the ice mixture is water we also account 
for the fraction of C2H4 or C3H6 ice compared to water ice 
in order to avoid overestimating our rate. We consider that 
not all the oxygen ice will react with the ethylene or propene 
ice because of the presence of other ice components in be- 
tween. We hence estimate a final value of 6xl0~ 15 cm 3 s _1 
at 20 K for the reaction between O and C2H4 and a rate 
of 5xl0 _14 cm 3 s _1 at 20 K for the reaction between O and 
C3H6 onto the grains. 

In order to reveal the contribution of these rates to the 
abundance of C-C2H4O and C-C3H6O in the ISM, we initially 
did not include any gas-phase pathway for the formation of 
these species. 

The percentage of depletion of oxygen atoms on the 
grains is certainly one of the critical parameters in the for- 
mation of ethylene and propylene oxide; Caux et al. (1999) 
found that the fractional abundance of oxygen (relative to 
hydrogen nuclei) on the grains is of the order of 10 -6 and 
this amount depends on the percentage of the oxygen that 
is depleted in the dense, cold cores. In our model, atomic 
oxygen freezes on the grains in the form of H2O, O and 
OH. It is known that oxygen containing species in interstel- 
lar ices consist mainly of H2O and also that the fractional 
abundances of water ice (relative to the total number of hy- 
drogen nuclei) are of the order of 10 -4 (Caux et al. 1999); we 
vary the percentage contribution of each channel by which 
oxygen containing species accumulate on the grains by vary- 
ing the efficiency of the different oxygen accretion pathways 
as shown in Table 2. 

Since we find 1 % of oxygen atoms in the ice (mO) is suf- 
ficient for the observational values of water ice to be satisfied 
we use this percentage subsequently for our model. Figure 1 
shows the fractional abundances of our selected species as a 
function of the time at the end of Phase I for a model with 
a final density of 10 7 cm" 3 . During this phase atoms and 
molecules freeze onto the grains and surface reactions can oc- 
cur due to the catalytic effects of the icy mantles. The dotted 
and dashed lines refer to the reactants, atomic oxygen and 
ethylene respectively, the solid line shows the ethylene oxide 
formation. Note that the fractional abundance of ethylene 
at the beginning of Phase I is low since it also mainly forms 
on the grains. We mostly focus on ethylene oxide because 



mO 



-13 - 
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Figure 1. Grain surface formation of selected species during 
Phase I. The fractional abundances (with respect to the total 
number of hydrogen nuclei) vary along with the time. 



results show that C-C3H6O is produced below the threshold 
of detection; this fact may mean that this molecule is not de- 
tectable in the ISM although we stress that we do not have 
a complete network of reactions (both in gas-phase and on 
grain surfaces) for c-CsHaO and its precursors. 

We conclude that the reaction between oxygen and 
ethylene on grains with a reaction rate of 6xl0" 15 cm 3 s _1 
(evaluated as described above) may be the main formation 
path for C-C2H4O. 

In order to compare our theoretical results with 
observations we ran further models where gas-phase re- 
actions were included despite the fact that the routes 
proposed to date (as mentioned in the Introduction) 
do not distinguish among the different isomers with 
formulas CbBUO. We then include in our model the gas- 
phase routes of formation for C-C2H4O listed in Table 
3. We insert the C 2 H40 species in the UMIST database 
with the reaction rates taken from the OSU database 
(http: / /www. physics. ohio-state.edu/~eric/research. html ) . 
Since the literature does not distinguish among the different 
isomers, we assume the ratios of 2:1.5:1 (as reported by 
Ikeda et al. 2001) for acetaldehyde, ethylene oxide and 
vinyl alcohol respectively to evaluate the rate of formation 
of each product. 

The different paths listed in Table 3 are mostly two- 
body reactions. Specifically we have charge exchange and 
ion-neutral reactions, electronic recombination of positive 
ions with electrons followed by the dissociation of the re- 
sulting molecule (dissociative attachment) and finally photo- 
processes where the dissociation or ionization of neutral 
species is performed by UV photons from a standard inter- 
stellar UV field. In the latter case, the reaction rates depend 
on the visual extinction of the source so we do not expect it 
to be efficient at the high visual extinction of hot cores and 
hot corinos. 

Before testing the sensitivity of the C-C2H4O abun- 
dances to the different physical conditions of potential as- 
tronomical sources, we investigate another important pa- 
rameter for the formation of ethylene oxide: the amount of 
C2H4 available on grain surfaces. Despite the fact that C2H4 
ice has not been detected to date, several experiments (Hi- 
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Table 3. Gas-phase paths of 0C2H4O formation and destruction. The rates are taken from the 
OSU database and they have been scaled by considering the ratio 2:1.5:1 as in Ikeda et al. (2001). 
R indicates a reactant, P a product. 
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to reproduce the relevant C-C2H4O observations. In a sim- 
ilar manner as described above for the case of oxygen, we 
vary the branching ratios for the different channels involved 
in the relevant hydrogenation pathways (see Table 4) in or- 
der to determine the minimum concentration of C2H4 on 
grains needed for the formation of ethylene oxide. As shown 
in Table 4, we investigated the variation of the branching 
of C2H2 hydrogenation to mC2H4 (m denotes a molecule 
in an icy mantle) between 1% and 40% (Table 4) and we 
obtain mC2H4 fractional abundances (relative to the total 
hydrogen nuclei) of 9.31xl0 -11 and 2.53xl0 -9 respectively. 
We find the lower abundance of 111C2H4 still generates vi- 
able quantities of C-C2H4O and so we fixed the branching 
of reaction 4 to 1% in all further studies. We then investi- 
gated the effect of the role of accretion of gas- phase C2H 4 
to form mC 2 H 4 , a process which competes with hydrogena- 
tion of C2H 4 to yield mC2H6 (reaction 1 and 6, Table 4). 
Here we find that a branching of only 1% in the direction of 
mC2H 4 still results in a c-C 2 H 4 abundance in agreement 
with observations. 

In summary, the above investigations of the hydrogena- 
tion pathways showed that a branching of only 1% towards 
both reactions 1 and 4 in Table 4 yields C-C2H4O abun- 
dances (3.9xl0~ 12 cm -2 ) comparable with observations and 
these branching ratios were used in all further modelling. 
Finally, we qualitatively model a sample of high-mass star- 
forming regions by looking at the influence of the physical 
parameters of the sources on the ethylene oxide abundances. 
Models are discussed in Section 3. 



Table 4. The main accretion and hydrogenation pathways for 



C2H2, C 2 H 4 , C 2 H 6 on 


grain surfaces. 


Reaction 


Branching ratios (%) 


1 C 2 H 4 -> mC 2 H 4 


10-1 


2 C 2 H 6 -> mC 2 H 6 


100 


3 C 2 H 2 -> mC 2 H 2 


10 


4 C2H2 mC 2 H 4 


40-1 


5 C 2 H 2 -» mC 2 H 6 


50-89 


6 C 2 H 4 -> mC 2 H 6 


90-99 



raoka et al. 2000; Kaiser & Roessler 1998; Barratta et al. 
2003; Bennett et al. 2006) have looked at its formation on 
grains under physical conditions relevant for the interstel- 
lar medium. Using different experimental techniques, C2H4 
was observed as a product following irradiation of CH4- 
containing icy mixtures and C2H4 was also predicted to be 
an intermediate in the hydrogenation of C2H2 (Kaiser & 
Roessler 1998; Bennett et al. 2006). In these cases, ethylene 
formation competes with the formation of the C2H6 that 
was found the most abundant molecule in the product mix- 
ture. In our model we include the hydrogenation of C2H2 on 
grains in Phase I to form C2H4 as well as C 2 H6 (see Table 
4). 

Based on the experimental studies performed to date 
(Barratta et al. 2006), each time C2H6 is produced, a small 
experimental signature that has been assigned to icy ethy- 
lene can also be observed; it is therefore appropriate to inves- 
tigate the minimum amount of ethylene on grains necessary 
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3 ASTRONOMICAL IMPLICATIONS 

Based on the study by Ikeda et al. (2001) ethylene oxide 
has been detected towards the Orion compact ridge and 
several hot cores. We qualitatively model these regions with 
UCL_CHEM by varying the final sizes and densities of Phase 
I and the temperatures in Phase II. In Phase I the non- 
thermal desorption is included as in Roberts et al. (2007). 
During Phase II icy mantles sublimate as in Viti et al. 
(2004) . In order to investigate the sensitivity of the chem- 
istry to the physical parameters of dense cores we run a grid 
of 16 chemical models by varying: (i) the size of the core, 
(ii) the final density of the collapsing core, (iii) the temper- 
ature, (iv) the adsorption energy (controlled by e) and (v) 
the percentage of the accreted species onto the grain sur- 
faces. Table 5 lists the different models and the C-C2H4O 
fractional abundances (/) (relative to hydrogen nuclei) ob- 
tained as outputs (see the last column). We also derive the 
column density (N) by using the formula below: 

N(c~C 2 H 4 0) = X x A v x N(H 2 ), (1) 

where X is the ethylene oxide fractional abundance, A v 
is the visual extinction and N(H2) is equal to 1.6xl0 21 cm -2 
which is the hydrogen column density in 1 magnitude. 



3.1 Sensitivity to variation in the density, 
temperature and size core 

Models from Ml to M10 reproduce the typical physical pa- 
rameters for an hot core model as described in Lerate et al. 
(2008); models M11-M16 best represent the characteristics 
for a compact ridge where we set a temperature of 80 K in 
order to look at the molecular abundances when the water 
desorption has not yet occurred. We consider a star of 25 
M© and 0.02-0.06 pc in size, characterized by a visual ex- 
tinction of about 300-500 mags. We then vary the density of 
the core (n H ) in the 10 6 -10 8 cm 3 range and we also modify 
the temperature (T) according to the observations towards 
the different astronomical sources, even if for T ^ 100 K 
we do not expect any influence of this parameter since, as 
already explained, the sublimation of ethylene oxide occurs 
simultaneously with the water desorption (at around 100 
K). 

If we compare, for instance, Ml, M5 and M6, no no- 
table differences in the column densities can be found: with 
increasing density the values are slightly higher as expected 
from the faster chemistry in a denser medium. On the other 
hand, keeping all the other parameters unchanged, the core 
size may influence the abundances of C-C2H4O: the final o 
C2H4O column density calculated from the output of M9 
is double the value estimated for Ml; the differences in the 
ethylene oxide abundances can vary up to one order of mag- 
nitude (see Ml and Mil). This observation is explained by 
the fact that the larger the core the smaller its visual ex- 
tinction (A v ) at a certain density. The magnitude of A v is 
related to the ability of the UV flux penetrating into the 
core and the UV flux leads in turn to the photodissocia- 
tion of species; in other words, in the larger regions we will 
observe a greater efficiency of molecular photodissociation. 



3.2 Sensitivity to variation in the percentage of 
the accreted species onto the grain surfaces 

Based on the study by Roberts et al. (2007) we introduce 
three non-thermal desorption mechanisms during the cold 
phase of our model (Phase I). Specifically, desorption result- 
ing from: H2 formation on grains, direct cosmic ray heating 
and cosmic ray-induced photodesorption. In particular, the 
surface formation of molecular hydrogen releases an energy 
of around 4.5 eV capable of inducing the thermal desorp- 
tion of molecules from the grains. This is a selective process 
which involves mainly the most volatile species. The cosmic- 
ray induced photodesorption is also a selective process that 
accounts for the number of molecules released to the gas- 
phase after each cosmic ray impact. These cosmic rays are 
able to ionize and excite the surrounding gas resulting in the 
emission of a UV photon, even where the visual extinction is 
high enough to extinguish any UV radiation propagating in 
from the outside of the cloud. This internal UV flux impinges 
on the grain surfaces and dissociates the absorbed molecules 
(mainly H2O). As for the case of H2 formation, molecules 
desorb due to energy released by recombination of OH and 
O; this type of desorption is non-selective. All the three non- 
thermal desorption processes are found to be very significant 
in dark cloud conditions where temperatures are low enough 
for the thermal desorption to be insignificant. Finally, direct 
cosmic rays are an efficient heating source able to penetrate 
deeply into molecular clouds; indeed when these energetic 
particles impact on the grains, molecules (adsorbed on the 
surfaces) start to desorb because of the energy transfer from 
the cosmic rays to the icy mantles. 

The efficiency of the non-thermal desorption of 
molecules depends on the percentage of the accreted species 
on the grain surfaces, which is in turn linked to the rate 
of freeze-out. The latter quantity is a free parameter in our 
models controlled by the product of the sticking coefficient 
and the average grain surface area. We vary the efficiency 
of freeze out (fr) in order to obtain different percentages of 
solid CO within a range of 75-98 % for the hot core models 
and 40-98 % for the compact ridge model. We finally in- 
vestigate the effects of varying the adsorption energy that 
is controlled by the e parameter. The latter is defined as 
the number of species removed each time an H2 molecule is 
formed. The value of e is uncertain and may be greater than 
1, but, as stated by Roberts et al. (2007), at large values, the 
desorption will be so efficient as to prevent the build-up of 
any mantle material until very late times. This justifies the 
choice of small values for e, 0.01 and 0.1, in order to look at 
the effect of this free parameter on the desorption rates. As 
reported in Table 5, the column densities are slightly higher 
when the H2 formation becomes the dominant non-thermal 
desorption mechanism (e = 0.1). 

In general, we find that all the changes made to the 
physical parameters (see Table 5) produce fluctuations of 
the estimated ethylene oxide column densities within two 
orders of magnitude; indeed our theoretical column densi- 
ties are in the 10 12 -10 14 cm -2 range for all models. This 
observation may seem surprising for the case of the com- 
pact ridge models since at 80 K (the temperature adopted 
for these models) the co-desorption of molecules with water 
has not yet occurred. This explanation for the constancy of 
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Table 5. List of chemical models and their parameters: size, density (ng), gas temperature (T) 
during Phase II, threshold adsorption energy (e), efficiency of the freeze-out (fr) during phase I, 
the percentage of mantle CO (mCO) given by the freeze-out parameter at the end of Phase I 
of the chemical model. The last column lists the C-C2H4O theoretical fractional abundances (/) 
(relative to the total number of hydrogen nuclei) obtained as outputs from each model. We also 
report its calculated column densities (N). 
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Figure 2. Fractional abundances of ethylene oxide in both gas- 
phase and grain surfaces as a function of time during Phase II. 

the column density of ethylene oxide may reflect the fact 
that non-thermal and thermal desorption provide the same 
contribution to the molecular abundances, but under differ- 
ent physical conditions. Specifically, the chemistry of denser, 
but smaller regions is mostly influenced by high temperature 
effects. However, in low density, but more extended sources, 
the chemistry is driven mainly by non-thermal desorption 
phenomena. 



4 RESULTS AND DISCUSSION 

Figure 2 reports the evolution of the fractional abundances 
of ethylene oxide in both gas-phase (solid line) and solid 
state (dashed line) during Phase II for Model M7 (see Table 
5). At the end of the collapse, C-C2H4O is mostly on the 
grain surfaces, but its abundance drops rapidly at around 
100 K when its co-desorption with water occurs (after 10 4 



Figure 3. Theoretical fractional abundances of three typical trac- 
ers of hot cores: CH3OH (solid line), HCOOCH3 (dashed line) and 
CH3CN (dotted line). 



yrs). At the same time we observe a small growth in the 
gas-phase C-C2H4O abundances thanks to the thermal des- 
orption of solid C2H4O. The amount of ethylene oxide in 
the gas-phase remains high at all times due to non-thermal 
desorption effects. 

We now compare the theoretical column densities with 
the observations towards several high-mass star-forming re- 
gions where C-C2H4O has been detected. In general, we 
achieve a good agreement, within an order of magnitude, 
with all of the observational values for C-C2H4O reported 
in Table 3 by Ikeda et al. (2001). Referring to the range of 
column densities reported by these authors, we consider M7 
as the model representing most of the physical characteris- 
tics of the regions listed in Ikeda et al. (2001). Results are 
reported in Table 6 where the calculated and observed col- 
umn densities for three typical tracers for hot cores: CH3OH, 
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Table 6. Comparison between our theoretical column densities (in cm 2 ) with those observed 
towards several high-mass star-forming regions. 



C-C2H4O CH3OH HCOOCH3 CH3CN 

Model 6.2xl0 13 2.4xl0 17 5.6xl0 15 2.2xl0 16 

Observations 1.3xl0 13 -1.7xl0 14a 1.4xl0 16 -1.4xl0 17a 1.0xl0 14 -1.2xl0 16a > 2.1xl0 16i > 



a Molecular column density taken from Ikeda et al. (2001) 
" Molecular column density taken from Wilner et al. (1994) 



HCOOCH3 and CH3CN are listed. The first two species were 
detected by Ikeda et al. (2001), the latter molecule was in- 
vestigated by Wilner et al (1994) which mapped the methyl 
cyanide lines in the millimeter spectrum of Orion KL region. 
The trends of the fractional abundances of these species as 
a function of time are displayed in Figure 3. The C-C2H4O, 
CH3OH (solid line), HCOOCH3 (dashed line) and CH3CN 
(dotted line) theoretical abundances perfectly match their 
estimated column densities. 

This result suggests that the reaction between atomic 
oxygen and ethylene on the grains is a suitable formation 
route for ethylene oxide. 



5 CONCLUSIONS 

We have investigated the feasibility of the reaction between 
oxygen ice and ethylene ice as possible route for the for- 
mation of ethylene oxide. For this purpose, we employed a 
chemical model for a collapsing core by updating the rate 
for the formation of ethylene oxide by using the recent ex- 
perimental determination of this quantity of Ward & Price 
(2011). We also tested the sensitivity of ethylene oxide frac- 
tional abundances to the changes in the main physical pa- 
rameters of the collapsing core. We finally reproduced the 
effects of non-thermal and thermal desorption of ethylene 
oxide during both the collapse and the warm-up phase re- 
spectively. Our work shows a good agreement between our 
estimated column densities for ethylene oxide (under hot 
core conditions) and those observed towards several high- 
mass star-forming regions (Ikeda et al. 2001). In order to 
verify the reliability of our comparison we also looked at the 
theoretical and observational abundances of other molecules, 
such as CH3OH, HCOOCH3 and CH3CN, known as good 
tracers of hot cores. We conclude that C-C2H4O can form 
on the grain surfaces thanks to the reaction between oxygen 
and ethylene. We stressed that only 1% of oxygen atoms de- 
pletion has to be in form of oxygen ice in order to produce 
ethylene oxide on the icy mantle (see Section 2); we also 
note that there is experimental evidence for a small quan- 
tity of ethylene being produced during the hydrogenation of 
acetylene to form ethane; therefore it is robust to assume 
that 1% of acetylene is hydrogenated to solid ethylene as an 
initial condition of our model. This amount of ethylene is 
already enough for the observed C-C2H4O column densities 
to be reproduced by our theoretical results. Further exper- 
imental studies are perhaps needed to extend our study to 
other cyclic species that have been detected in the ISM, al- 
though the complexity of larger molecules may prevent their 



isomeric speciation as is the case for the recent laboratory 
study of the formation of CsHgO. 
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